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abstract 

The aim of this paper is to compute transitions amplitudes in quantum cosmology, 
and in particular pair creation amplitudes and radiative transitions. To this end, we 
apply a double adiabatic development to the solutions of the Wheeler-DeWitt equation 
restricted to mini-superspace wherein gravity is described by the scale factor a. The 
first development consists in working with instantaneous eigenstates, in a, of the matter 
Hamiltonian. The second development is applied to the gravitational part of the wave 
function and generalizes the usual WKB approximation. We then obtain an exact equa¬ 
tion which replaces the Wheeler-DeWitt equation and determines the evolution, i.e. the 
dependence in a, of the coefficients of this double expansion. When working in the grav¬ 
itational adiabatic approximation, the simplified equation delivers the unitary evolution 
of transition amplitudes occurring among instantaneous eigenstates. Upon abandon¬ 
ing this approximation, one Ends that there is an additional coupling among matter 
states living in expanding and contracting universes. Moreover one has to face also the 
Klein paradox, i.e. the generation of backward waves from an initially forward wave. 
The interpretation and the consequences of these unusual features are only sketched in 
the present paper. Finally, the examples of pair creation and radiative transitions are 
analyzed in detail to establish when and how the above mentioned unitary evolution 
coincides with the Schrodinger evolution. 
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1 Introduction 


In General Relativity, one loses the usual Hamiltonian structure in which time parame¬ 
terizes the evolution since the theory is invariant under arbitrary changes of coordinates. 
In order to understand the consequences of this loss in quantum cosmology, one must 
begin by reexamining the status of time in quantum mechanics and quantum held theory. 
In particular, we draw the attention to the distinction between the adiabatic deforma¬ 
tion of the states and non-adiabatic transitions amongst states. This distinction will be 
very useful upon addressing the question of the evolution in quantum cosmology. 

For simplicity, we first consider the case in which the Hamiltonian is time inde¬ 
pendent. Then, its eigenstates are stationary and, in the absence of interactions, the 
probability amplitudes to find the system in each eigenstate are constant. When a time 
independent interaction Hamiltonian is added, its sole role is to displace the energy levels 
from their values defined by the free Hamiltonian to their dressed values defined by the 
(stationary) eigenstates of the full Hamiltonian. 

When the Hamiltonian is time dependent, it has a double role. First, it modifies 
adiabatically the energy levels defined by its instantaneous eigenstates. In the limit of 
slowly varying situations, one recovers the former case in which the amplitudes stay 
constant. Secondly, it induces non adiabatic transitions between instantaneous eigen¬ 
states, a truly time dependent phenomenon since it characterizes quantum jumps. The 
importance of these transitions is controlled by d t AE / AE 2 where A E{t) is the energy 
difference between the instantaneous eigenstates. It goes without saying that as long as 
the interaction Hamiltonian is hcrmitian, the existence of these transitions do not re¬ 
quire a revision of the interpretation of quantum mechanics since the norm of the state 
stays constant. 

In quantum held theory, time dependent gravitational fields induce both effects, 
although they are somewhat hidden when one works -as is usually done- in a Heisenberg 
picture. For instance for a free scalar held one would proceed by computing the solutions 
of the Klein Gordon equation. Its WKB solutions describe the adiabatic deformation of 
positive frequency waves and lead to no transitions in Fock space. Instead, backscattering 
between positive and negative frequency waves correspond to non-adiabatic transitions 
between states, i.e. pair creation from vacuum. What plays the role of A E(t) in quantum 
mechanics is 2u (f), the energy of a pair. Note that the existence of the backscattering 
amplitude leads to the “Klein paradox”, namely the increase of the current carried 
by the forward wave when some backward wave is generated (since the total current 
is conserved). This imposes that one abandon a hrst quantized framework, since no 
conserved positive local quantity can be built from these solutions. 

In the hrst part of this article we review the notions of adiabatic evolution and 
adiabatic transitions in a cosmological context in which the evolution of the scale facor 
a(t ) is given from the outset. We show that a itself can replace time and be used as the 
parameter which controls both the adiabatic deformations of the states and their non 
adiabatic transitions, although time still appears in the relative phase accumulated by 
neighboring states. This prepares for the analysis in quantum cosmology where a is a 
dynamical variable. 
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We illustrate this formalism by two examples. The first is the recovery of the Golden 
Rule, i.e. that transition probabilities grow linearly with proper time lapses. We show 
how the residual time dependence mentioned above guarantees that the transition rate 
is constant. (Thus a collection of such systems can be viewed as a simple realization of a 
“clock”.) The second example, taken from field theory, is the process of pair creation in 
cosmology. We show that pair creation acts are described, in the Schrodinger picture, by 
non adiabatic transitions between states having different instantaneous particle number. 
Moreover, this second example introduces uscfull tools to handle the gravitational part 
of the solutions of the Wheeler-DeWitt equation. 

In the second part this article, we work in quantum cosmology. In that framework, the 
physical states are solutions of the Wheeler-DeWitt (WDW) equation. Indeed, instead 
of having as before a Hamiltonian framework in which processes are parameterized in 
time, one has a constraint that expresses the invariance of the theory under arbitrary 
changes of time (of time only since our analysis is restricted to mini-superspace) [|lj]. 
Therefore, to make contact with the conventional analysis, one must understand how 
to define adiabatic deformations of matter states and how to compute non adiabatic 
transition amplitudes without having a physical time parameter at our disposal. 

The first point is easily addressed. Indeed, since the Wheeler-De Witt equation is 
the sum of a kinetic term for gravity, a potential term coming from the curvature of 
space, and the matter Hamiltonian, Hm(cl), the adiabatic matter states are simply the 
instantaneous -in a- eigenstates H M (a). Therefore they coincide with the instantaneous 
eigenstates -in t- defined in the Schrodinger framework, since in cosmology, the whole 
dependence in t of H M comes through a(t) only. 

To address the second point, it is necessary to carry out a second adiabatic expansion 
applied to the gravitational part of the wave function. This generalizes the usual WKB 
expansion which has been used in [@ 0 @ @ to extract a time parameter which governs 
the Schrodinger equation in the mean geometry. The central result of this paper is the 
(exact) equation satisfied by the coefficients of this double adiabatic expansion for both 
matter and gravity. 

When one works in the adiabatic approximation for gravity, this equation has a 
very simple structure: it is a linear, first order equation in d a and bears a simple in¬ 
terpretation. Indeed, the identification of what replaces the probability amplitudes in 
Schrodinger equation is straightforward. It describes therefore the non adiabatic transi¬ 
tions amplitudes of matter parametrised by a rather than t. In this, we generalize 0-0 
which were based on time dependent perturbation theory. Indeed, in the adiabatic 
treatment, the dressing energy of matter states is automatically included in the wave 
functions describing gravity and non-adiabatic transitions are evaluated from gravita¬ 
tional wave functions wherein this back-reaction effect is fully taken into account. Thus 
we no longer meet the problem 0 0 of defining transitions amplitudes that are non lin¬ 
ear in the coupling constant. These advantages are illustrated in the last part of this 
article when we return to the examples of radiative transitions and particle creation 
now treated in quantum cosmology. We explicitly show when and how one recovers the 
background field result. Mathematically this requires to legitimize a first order expan¬ 
sion in the energy difference AE(a) that characterizes the transition upon investigation. 
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Physically it means that the universe be macroscopic 0 [|]], i.e. that the matter sources 
driving gravity be macroscopic. 

When one abandons the adiabatic treatement for gravity, one meets new effects which 
have no counterpart in quantum mechanics. These arise because one must consider 
superpositions of solutions describing contracting and expanding universes. In this one 
inevitably encounters the “Klein paradox”, i.e. the generation of backscattered solutions 
describing contracting universes when the wave was initialy purely forward. As in QFT, 
these backscattered solutions lead to an increase of the current carried by the forward 
wave describing expanding universes. Our dynamical equation takes exactly into account 
the contribution of these backscattered solutions. It can thus be exploited to analyze 
both conceptually and quantitatively their consequences. One must face indeed a double 
problem. First, one must adopt a framework in which the generation of backscattered 
waves can be interpreted. In this respect, one must probably proceed to a “third” 
quantization. Secondly, there is also the more pragmatic problem which consists in 
determining their consequences onto the matter propagation restricted to the expanding 
sector, i.e. onto matter belonging to “our” universe. 

^Frorn our exact dynamical equation, the much debated question of the unitarity of 
the evolution can now be phrased in precise mathematical terms. One should indeed 
first identify what mathematically corresponds to the probability amplitude to find a 
given adiabatic state at a certain a. Only then can one hope to settle the question of 
the unitarity of their evolution. In this article we shall only sketch these problems. We 
hope to return to them in a future publication. 

We conclude by noting that the analysis of the Wheeler-De Witt equation presented 
here possesses illuminating analogues in other physical systems. Indeed the basic prob¬ 
lem, which is not specific to quantum cosmology, is to enlarge the dynamics so that a 
degree of freedom which was treated classicaly (the radius of the universe a(t )) is treated 
quantum mechanicaly. Some of the conceptual issues raised by this enlargement will of 
course be specific to each problem. 

The first analogue is a uniformly accelerated detector [p] which can make non adi¬ 
abatic transitions from its ground state to its excited state, and vice versa. (The de¬ 
scription of these events as non adiabatic transitions can be found fllQ|] . The treatement 
presented there is also closely related to our analysis of the Golden Rule). These transi¬ 
tions are characterized by a temperature, in similar fashion that the pair creation events 
in cosmology which we study in this article are characterized by a temperature. In these 
works the accelerated system follows a given classical trajectory. However, in order to 
take into account recoil effects one must go beyond this background field approxima¬ 
tion, and describe the accelerated system by WKB wave functions|ll]] or even by second 
quantized fields[|T2|. Then not only is the trajectory treated quantum mechanically, but 
one can also produce pairs of accelerated systems, which is strongly reminiscent of third 
quantisation in cosmology. In both cases what should be interpretated as the proba¬ 
bility amplitudes to find a given adiabatic state in the forward sector do not obey a 
closed equation. In other words, “our” expanding universe is not isolated. The accel¬ 
erated detector system also exhibits a remarkable thermodynamic consistency between 
the radiative transitions of the detector and the pair production events. In quantum 
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cosmology we will also find an intimate connection between the dynamics of transitions 
in the universe and the backscattering of universes. We have however only exhibited 
this connection in this paper, and have not investigated in detail its consequences. 

The second analogue is the description of electronic non elasticity, i.e. transitions 
between electronic states, during atomic or molecular collisions |]~5| ||14||. The simplest 
treatment is a background field approximation wherein the nuclear trajectories are given 
classical functions R(t). The electrons then evolve according to the time dependent 
Hamiltonian H e i(R(t)). (This corresponds of course, in our case, to QFT in a given 
space time.) However one can go beyond this approximation and treat the nuclei quan¬ 
tum mechanically (although in first quantization) ||T5| |16| . To this end one carries out a 
double adiabatic development in the electronic adiabatic states and the nuclear WKB 
wave functions. The equation for the coefficients of this double expansion describe the 
non adiabatic transitions of the electrons, parametrised by the nuclear coordinates in 
exactly the same way that the equation for non adiabatic transitions of matter is re¬ 
covered in quantum cosmology. This equation can be used to investigate the validity of 
the background field approximation. In particular one can investigate in this context to 
what extent the notion of a (temporal) succession of electronic transitions depends on 
the nuclei being in a tight wave packet. The analogy with quantum cosmology obtains 
therefore both at the formal and at the conceptual level (save for the interpretation of 
the backscattered waves). Hopefully these systems could be used to investigate experi¬ 
mentally some of the conceptual problems of quantum cosmology. 

As a final comment, we wish to acknowledge the fact that many people already ap¬ 
peal to the use of adiabaticity to investigate matter dynamics in quantum cosmology. 
Suprisingly, to our knowledge the full power of the adiabatic treatement (by opposition 
to the adiabatic approximation ) has not been exploited previously. However, after con¬ 
cluding this work, we became aware of refs. 00 wherein the full treatement is used, 
although the motivations of these authors and the chosen exemples differ significantly 
from those of the present paper. Indeed they have restricted their analysis to the “Hawk¬ 
ing model” wherein matter consists only of a massless scalar field (because of its possible 
relevance to inflation), whereas we insist on its generality and applicability to all matter 
processes. Furthermore they have not used this formalism to analyse the origin of time 
in quantum cosmology and the related conceptual issues which is the main focus of the 
present paper. 


2 The adiabatic treatment in quantum mechanics 

2.1 Instantaneous eigenstates and non adiabatic transitions 

In this section we recall the basic features of the adiabatic approximation applied to 
quantum mechanics. We also review how to compute non-adiabatic transitions. To 
illustrate the usefulness of this formalism in a cosmological setting, we apply it in sub¬ 
sections 2.2 and 2.3 to matter transitions and pair creation amplitudes. We put special 
emphasis in the role played by the expansion of the universe described by a = a(t ) 
since, in cosmology, the expansion is the cause of non adiabatic transitions. Indeed, in 
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the absence of expansion, the matter would be described by stationary (un-interfering) 
eigenstates. This sets the stage for Section 4 where the same transitions will be described 
in quantum cosmology, wherein gravity is described by wave functions, solutions of the 
Wheeler-De Witt equation. 

The Schrodinger equation that governs the evolution of the state |is 

id t \jp(t)) = H(t)\jp(t)) (1) 

We expand \i^(t)) in terms of the instantaneous eigenstates of the time dependent Hamil¬ 
tonian H(t): 

IV’W) ='52c n (t)\‘tl>n(t)) (2) 

n 

where 

H(t)\il> n (t)) = E n (t)\ip n (t)) 

(^ n (t)|^ m (f)) = 8 nt m (3) 

Since H{t ) is hcrmitian, the eigen-energies E n (t) are real. Inserting the decomposition 
eq. (H) into eq. (H) yields the equation for the c n (t ) 

id t C n (t) = E n (t) C n (t ) + i J2 (dti>m(t)\'tpn(t)) C m (t ) = H nm (t) C m (t) (4) 

m m 

where we have introduced, for later convenience, the effective Hamiltonian H e ^(t). Note 
that the antisymmetry of ( dtif) m (t)\if) n (t )) = and the reality of E n (t) 

ensure that H e ^(t) is hcrmitian. Therefore the conservation of probability J2 n |(t) | 2 = 
1 is guaranteed. 

Since the states \ij; n (t)) are defined up to an arbitrary phase, we may use this freedom 
to require 

(' l l> n (t)\d t ip n (t)) = 0 (5) 

Note however that if \^ n {t)) evolves cyclically during the interval f 0 , then the state at 
t\ will in general differ from state at time to by a phase, known as the Berry phase fT9|. 
Berry’s phase will not be of concern to us in the present paper since in mini-superspace 
the matter states depend only on one parameter, the scale factor a. In this case Berry’s 
phase vanishes, since all paths can be mapped one onto another by a reparameterisation. 
However if one goes beyond mini super space, one should take the Berry phase into 
account. 

Upon absorbing, as in time dependent perturbation theory, the diagonal part of 
H e lf (t) into the definition of new coefficients c n (t) = c n (t)e l J dtEn ^\ one obtains 

d t c n = £ Cm(t) (6) 

m^n 

At this point, it is important to notice that this equation is “almost” reparametrisation 
invariant. The only term that does not satisfy this invariance is the phase containing the 
change in energy E m (t') — E n {t'). As we shall show in section 2.2, it is this factor that 
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guarantees that the Golden Rule is recovered, i.e. that probabilities increase linearly 
with respect to (proper) time. In a cosmological situation, the sole time dependence of 
H{t) comes through the expansion law a = a(t). Therefore the eigenstates \ijj n (a)) and 
the eigenvalues E n (a ) depend on t through a(t ) only, by virtue of their definition, see 
eq. (|3|). Thus, one can write eq. (|B|) directly in terms of a as 

daCn = E {daMaMM) e _< (W(iW-W) Cm(a) (7) 

rriy^n 

We shall see in Section 3 that it is in that guise, i.e. wherein t appears only as a factor 
of E m (a) — E n (a) and in the form dt/da, that the evolution of the c n will be delivered 
in quantum cosmology. 

To proceed we must evaluate ('i/} n (t)\d t 'il) m (t)) ( n j- m). If H is non degenerate this 
is given directly in terms of the Hamiltonian. Indeed differentiating eq. @ 

dtH\il> n ) + H\dti/> n ) = d t E n \i/i n ) + E n (t)\dM ( 8 ) 


and taking the scalar product with (^ m | yields 




(V> m\dtH\i> n ) 


E n ~ E 


m 


n 7 ^ m 


( 9 ) 


Thus the time dependence of H(t ) leads to non-vanishing off diagonal elements in Ed e E. 
These have two physical consequences. First they modify the instantaneous eigenstates 
and eigenvalues |20|. |2T| . Secondly, they induce transitions among eigenstates 

The first effect is described by carrying out an instantaneous diagonalisation of 
H e E(t) to obtain new eigenstates and eigenvalues: H e E[t]) = E^(t)^^(t)). 
After having diagonalised iF e ^(t), one can decompose | ijj) as | ^(t)) = l^nKt)). 

The evolution of the new coefficients is given by a new effective Hamiltonian 

H e ^^\t). One can then iterate the procedure and diagonalise H e ^^ l \t). This gen¬ 
erates an asymptotic series for the state |The series starts diverging when the 
difference between successive terms becomes comparable to the amplitude for non adia¬ 
batic transition. 

To work in the adiabatic approximation consists in neglecting the second effect, 
namely non adiabatic transitions between instantaneous eigenstates. In general, to com¬ 
pute the amplitude of these transitions, it is sufficient to consider only a couple of states, 
say |-0i) and l^)- Then eq. ([|) reduces to 

d tCl = (d t HMe-*!**&-*> ci 

d t c 2 = {d t ^ 2 ) Cl (10) 



Let us further assume that at t — —oo the system is in the state - 01 ; be. C\(—oo) = 1, 
c 2 (— oo) = 0. Then the transition amplitude is given by C 2 (+oo). To evaluate it, we 
proceed exactly like in perturbation theory by taking c\ = 1 in eq. O to obtain 

c 2 (+oo) ~ f dt(-02|cV0i) e~ z X-°° dt ( El( * ) _E2( *) ) ( 11 ) 

■' —OO 
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We consider only the case in which the resonance condition Ei(t) — E 2 (t) = 0 has no 
real solution. (This is generally the case in physical problems, because if their is an 
interaction between the two energy levels the degeneracy is lifted to give rise to an 


avoided crossing [p2| ). Then, when the transition amplitude is small, it can be correctly 
evaluated by a saddle point treatement, exactly like in other cases of classically forbidden 
transitions such as tunneling processes. The saddle is located at the comp lex value of t*, 
solution of E] (t*) — E 2 (t*) =0. The result is thus 


. * 

( - \ —i f dt'(Ei — E2) 

C2(+oo) 


( 12 ) 


A rigorous analysis shows that C tends to 1 in the adiabatic limit, d t (Ei — E 2 )/(Ei — 

TP ~\2 


0, see [|23|. If the problem cannot be reduced to that of two coupled states, 


then the analysis is more complicated because the non adiabatic transitions can occur 
by “hopping” from one energy eigenstate to the other, see p4 . 


2.2 The Golden Rule and non adiabatic transitions 


In this section we apply the adiabatic formalism to radiative transitions in cosmology 
in order to display the crucial role played by the classical expansion a = a(t). The 
particular problem we consider is the disintegration of a heavy state of mass M into a 
light state of mass m and a conformally coupled scalar quantum (f>. In the language of 
Unruh||, M and m represent the excited and ground state of a particle detector of the 
(j) field. As in the Unruh analysis, the recoil of the heavy system is neglected. Because 
this transition takes place in an expanding universe, the energy of the emitted quantum 
(measured in proper time) is k/a(t ) where k 2 is the conserved eigenvalue of the three 
Laplacian. The free initial and final energies are thus E® nitial = M and Ej inal = m + k/a. 
They coincide when the Doppler red shifted energy of the photon is equal to M — m, 
i.e. when 

k 

a{t 0 ) =a 0 = — (13) 

Thus, in time dependent perturbation theory based on free states, the transition will 
occur around that value of a wherein the energy of the photon resonates with AM[[7|. 
In what follows, we shall see that non adiabatic transition will also occur in the vicinity 
of to even though the meaning of “to make a transition” will change. 

The unperturbed states are coupled by an interaction Hamiltonian 


H int (a ) = g (| M >< m\4> + Ac.) = g [\M >< m\ ^ ^^ (b k + b[) + h.c. 


Thus the full Hamiltonian is 

H(a ) 


( M g/ak 1 ' 2 \ 
\ g/ak 1 / 2 m + k/a J 


(14) 


(15) 


As pointed after eq. the time dependence comes only through a = a(t). Therefore, 
one can directly work with a as relevant parameter. Because of the non-diagonal elements 













induced by H int , the degeneracy of the free states at a 0 is lifted to give an “avoided 
crossing” 11 when one considers the instantaneous eigenstates, | i[} + ) and |-0_), which 
diagonalize H(a). Indeed, their corresponding eigenvalues are 


E±(a) }l ' ' k/a ± \[{M -m- k/a) 2 /A + g 2 /a 2 k (16) 

They coincide only for the complex values 

a ± = a ° ± A mIi/2 = ± iAa i ( 17 ) 

The principle effect of this avoided crossing (for real values of a) is that the instantaneous 
eigenstate \if; + ) (|V 1 -)) which for a -C ao coincides with the free initial (final) state 1 4’initai) 
(I ^final)) tends to the free state \ipfinal) (| ^initial)) for a > a 0 . Thus what corresponds 
to the transition amplitude in terms of free states is given by the amplitude to stay in 
the same instantaneous eigenstate, whereas the amplitude to stay in the same free state 
is given by the amplitude to make a non-adiabatic transition from |^ + ) to |-0_). 

By writing the probability to make a non-adiabatic transition, that is to stay in the 
excited state M as 


P 


no transition 


(k) = 


- c - 2 Ik 


(18) 


using eqs. (0, 0), we have, 

Tk = Im 


r a X dt 

/ da —(a) [E + (a) — E_(a)} 
ao da 


(19) 


since only the integral for imaginary (positive) values of a contributes to Ik- To evaluate 
this integral, we assume that g be small enough so as to justify the following approxi¬ 
mation 


Ik — 


1 dt 
ao da 


( a 0 )Im f da\J (A Ma — k) 2 + 4 g 2 /k = 
J ao 


ng 2 dt 
k 2 da 


(ao) 


( 20 ) 


Within this quadratic approximation, the probability for no transition is given exactly 
by eq. (|TSp, see [|nj]. This exact result plays an important role here since it guarantees 
that the probability not to decay decreases exponentially in (proper) time, see eq. 

This probability is given by 


Pno transition (fo 11 Pno transition ( k ) e 

k(a) 


- 2 / 
- tD J 


k ( a ) d 3 k 


(2 ttY 


Ik 


( 21 ) 


To carry out the integral over k in the exponent, we first use the isotropy of the space-time 
to write d 3 k = 4:7ik 2 dk (for simplicity, we work with flat three-geometries). Furthermore, 
by using the resonance condition eq. O. we can write dk as AMdao and thus relate 
the final value of a to a maximum value of k , see for a more detailed treatment 
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concerning the relationships between the integration over momenta and an integration 
over time (or a). Using that relation and eq. (|0]), one gets 


p 


u 2 am r a 


( 22 ) 


corresponding to an exponential decay probability in proper time. It agrees with the 
probability evaluated using more traditional techniques^. 

It is interesting that in the present formalism the role of time and energy seem to be 
reversed compared to usual treatment of the Golden Rule. Namely if the particle does 
not decay it seems to be propagating in complex time, see eq. (0)> whereas one usually 
ascribes a complex energy to unstable states. And the time T = f dt which multiplies 
the decay rate in eq. (|2^) arises from an integral over fc, whereas it is usually obtained 
by letting the excited states interact with the field through Hi nt for a finite time T. 

The treatment presented here can be enlarged to include the effects of momentum 
recoils. This gives rise, as discussed in [7J, to the concept of spatial displacement. 


2.3 Particle creation as a non adiabatic process 

In this section we apply the adiabatic formalism to particle creation in the universe. 
This particle production |26|, [27] is a generic feature of quantum field theory which is due 
to the expansion law a = a{t). 

When the cosmological expansion is described by a classical law, the equation for 
free matter fields is linear in the dynamical variables. Therefore, it is particularly useful 
to work in the Heisenberg picture since probability amplitudes to create pairs can be 
directly evaluated from the sole knowledge of the Bogoliubov coefficients. However, 
upon abandoning the background field approximation for gravity by working in quantum 
cosmology, one looses the linearity of the field equation since a is now a dynamical 
variable coupled to matter at a quantum level. As noted in |(|. this loss of linearity 
imposes to work in a Schrodinger picture since every transition amplitude should be 
computed separately. Therefore, in what follows, as a preparation for the analysis in 
quantum cosmology, we study pair creation amplitudes in the Schrodinger picture. In 
this picture, they correspond to non-adiabatic transitions among instantaneous eigen¬ 
states. 

We shall work in a closed Robertson-Walker universe 


ds 2 = —dt 2 + a 2 (t) ^2 hijdx l dxi 


where 


dr 2 




hij = -- ~ + r 2 d9 2 + r 2 sin 2 9d(j) 2 


13 1 _ r 2 

The action for a massive hermitian field is 


(23) 


(24) 


L = - J dt J d 3 xVha 3 


<9 f 0<9 t 0 — a 2 h l:, di(j)dj(l) — m 2 (j) 2 


(25) 
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The symmetry of the space allows <f> to be decomposed as 


H(a) = J2 [«“ 3 n fc n_ fc + a 3 nl(a)M_ k \ = J2 H k {a) 

k k 


0(f, x) = y k (x)Mt) (26) 

k 

where 34 are normalized eigenfunctions of the spatial Laplacian. Using eq. (^SJ), the 
action reads 

L = [dta 3 (t) [d t (f) k d t (j)- k - (27) 

k J 

where kl k {a) = (m 2 + fc 2 a -2 ) 1 / 2 is the time dependent frequency at fixed k. 

The momentum conjugate to <p k is = a 3 dt4>- k , and the Hamiltonian is 


(28) 


We now carry out the instantaneous diagonalisation of H k (a). To this end we ex¬ 
press the Schrodinger operators (f> k and in terms of the a-dependent creation and 
destruction operators b k (a ) 

(b k = —=-— (bk + tfu) 

^/2(a 3 ^) 1/2 

n fc = -j=i{a 3 n k )V 2 {b- k -bl) (29) 

Inserting eq. (|29|) into eq. (|8|) yields 

H k (a) = tl k (a){b\b k + b ] _ k b_ k ) (30) 


where we have subtracted the ground state energy by normal ordering with respect to 
the instantaneous operators b k . The instantaneous eigenstates are therefore those of 
definite particle number 

(tf ) n k 

I Vw(«)> = n^lOa) 

k Vn k '. 

E {nk }(a ) = J2^k(a)n k (31) 

k 

where {n k } denotes the collection of occupation numbers and |0 a ) the instantaneous 
vacuum state. 

The equation governing non adiabatic transitions is 

d aC{n k } = E (3a4W}IVw) e ~^ ( ^ C{n k ,} (32) 

i n k'} 

where the state of the (j) field has been decomposed as 

l^(«)) = 51 c {n k }(a) 14w(a)) ( 33 ) 

{n k } 
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To obtain an expression for {'^{n k }\d a '^{n k ,})we must calculate d a H k , see eq. ([]), 


d a H k = dA (b\b k + bl k b_ k ) + n fc (a) ^3^ fc) (b\b ] _ k + b k b. k ) (34) 


d a H k contains two terms. The first is diagonal in the basis of instantaneous eigenstates 
and is due to the adiabatic deformation of the eigen-energy through kl k (a). The second 
term is due to the time dependence of the operators b k (a), b\(a). It couples states 
which differ by a pair of particles of opposite momentum and engenders non adiabatic 
transitions. Thus the corresponding physical processes are pair creation (or destruction) 
of particles with opposite momentum. 

We now consider the process of creation of a pair from vacuum. Thus we seek the 
amplitude for the non adiabatic transition from the ground state |0 o ) to the two particle 
state bltf_ k \O a ) = |lfc, l- k )- We assume that the creation process be small enough so as 
to legitimate a description restricted to these two states. In that case, the equations are, 
see eq. (0, 


d a CO 


daCk,—k 


1 9 „(o 3 n t ) jj 

2 a 3n t e 

1 d a { 

2 a?Sl t 6 C ° 


(35) 


Thus, for rare creation acts, the probability to produce a pair is, see eqs. (|12|), 
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pair | 


12 ~ Im 

4 [ tl ' a 1 y/m 2 +k 2 /a 2 dt 

= e 



(36) 


where a* is the complex saddle point defined by Vt k (a*) = 0 in the positive imaginary 
half plane. 

To illustrate the above formalism, we apply it to de Sitter space. In that case, the 
expansion is given by 

a(t) = —j=coshy/h.t (37) 

Thus the complex saddle point of eq. (|36l) is 


= ±ik/m 


t± = to ± iAti = —j=arcsinh — -j= ± i(/r/2v / A 


y/A 


kV A 


(38) 


It lies in the vicinity of the real time t 0 at which k/a ~ m, i.e. where the particle ceases 
to be relativistic. When m > v^A, the probability to produce a pair is small and eq. 
( |36| ) can be used. The integral is most easily evaluated in the limit of large k, i.e. for 
|a*| 1/m. Indeed, in that limit, as in eq. (|20|), the integrant is quadratic in a 


Im 


4 f Jm? + k 2 /a 2 dt 

<7 — 00 


= Im 


~ Im 


r da ^ dtlda) Vm 2 a 2 + k 2 

lao CL 


alV A ■ 


day/m 2 a 2 + k 2 


ao 


= 7 rm /VX (39) 
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where we have used da/dt = ay /A valid for large a, and replaced a 0 = a(t 0 ) = k/m. 
Thus we have _ 

\A pair \ 2 ~ (40) 

This corresponds to the Boltzmannian tail of a thermal distribution at the de Sitter 
temperature T = yfX/'K. This is as it should be since the use of the saddle point 
calculation and neglection of multi pair states cannot reproduce the full Bose Einstein 
distribution [[27] . 


Before considering how to recover these transitions in quantum cosmology, we shall 
first relate the present analysis to the more conventional analysis of pair creation based 
on the solutions of the Klein Gordon equation. The main motivation for this study is 
that it introduces tools that will be extremely usefull upon studying the solutions of the 
WDW equation. 


2.4 The WKB approximation and pair creation 

The aim of this section is to phrase in adiabatic terms the standard results concern¬ 
ing the WKB approximation and the backscattering amplitude above a barrier. This 
will clarify the relationship between pair creation probabilities as calculated above and 
backscattering. The necessity of normalizing properly the Wronskian of the WKB so¬ 
lutions will be stressed. This will also prepare for the analysis of the solutions of the 
WDW equation for which a correct normalisation will also be essential. 

Consider a one dimensional harmonic oscillator whose frequency depends on time 

<9 2 +u; 2 (f) (p(t) — 0 

The WKB solution of unit Wronskian is 

-if* 

P 


(41) 


dtu(t ) 


x(t) = 


2 u(t) 


One can then expand the exact solution in terms of y(t) and xif)* as 

<p{t) = c{t)x{t) + d(t)x*(t) 

We determine completely the coefficients c and d by requiring that 

= w(t) (c(f)x(f) - d(t)x‘(t)] 


(42) 


(43) 


(44) 


This guarantees that c and d are constant in the WKB limit, i.e. in the adiabatic limit 
d t uj(t)/uj(t) 2 —> 0. Moreover, the conserved current (or Wronskian) takes the simple 
form 

( p*id t (p = |c(f)| 2 — \d(t)\ 2 = constant (45) 

This equation replaces the unitary relation governing the coefficients c n (t), solutions of 
the Schrodinger equation, see 2.1. To obtain c and d, one differentiates eq. (fhi]) and 


compares the result with eq. ([44|) , so as to obtain 

Xd t c + X*dtd 


\ — (xc + X*d) = 0 

Z CO 


(46) 
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Now differentiating eq. ([44]), inserting the result into eq. ([U]), and using eq. ( f46|) to 
eliminate either d t c or d t d yields the coupled first order equations 


o 1 dfUJ —2i f 4 dtuj i 

o t c = -e a 

2 u 

8 t d = e +2i f dt “c (47) 

2 u 


These equations are equivalent to the original equation for ip, eq. (|4l|) . This has been 
discussed in detail in [28]. They constitute a convenient starting point for the evaluation 


of the non-adiabatic transitions, which in the present case correspond to backscattering. 
As in Section 2.1 they are given by the value of value of d at late times, given that 
c = l,d = 0 at early times. Using the same type of analysis as in section 2.1, one 
obtains 

,/ i \ f + °° A+^^ — 2i f* (j(t')dt' 

a(+oo) ~ / dt -e 

J —oo 2u 


(48) 


As before, the saddle point is at the complex value t* such that u>(t*) = 0. The norm of 
d is thus 

|d(+oo)| 2 ~ C e- 4Im ^-oo^ dt ' (49) 


A rigorous analysis shows that C tends to 1 in the adiabatic limit [ fZDj| . 

To obtain the relation between this norm and the probability to produce a pair of 
quanta as calculated in section 2.3, we must first determine the frequency u(t) that corre¬ 
sponds to Ufc(a) which caracterizes the instantaneous energy of a particle of momentum 
k. A standart calculation yields 


oj(t) = &, k (a(t)) 


3 (d t a) 2 3 <9 2 a\ 

4 a 2 2 a ) 


(50) 


The additional terms arise from the factor of a 3 / 2 which relates the amplitude of 4>k(t) 
appearing in eq. (|28|) to (pit). It is important to note that the backscattering amplitude 
and the pair creation amplitude are related by |d(+oo) | 2 = | A pair | 2 only if the Wronskians 
of the WKB waves that multiply c and d are chosen in conformity with the second 
quantized framework, i.e. are equal to one. In the present case of a single linear equation, 
only the relative normalisation plays a role, however upon introducing coupling among 
many oscillators, the normalisation is completely fixed, see [0,0-10,0. 

We wish to emphasize this last point since it will play a crucial role in quantum 
cosmology. The question of the choice of the Wronskians can be considered from two 
different perspectives. The usual approach is based on the second quantized framework. 
In that case, the choice of unit Wronskian is directly dictated by the particle interpre¬ 
tation of the Heisenberg creation and destruction operators. However, there is another 
point of view which is more related to our present understanding of quantum gravity. It 
is based on the semi-classical behaviour of some dynamical variable (say a). In physical 
circumstances in which a behaves “almost” classically, two different schemes can be used 
to describe quantum processes occuring to the other variables. First, one can work in the 
background field approximation for a (in which case, one has quantum fields in a given 
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gravitational background). Secondly, we can work in an extended quantum framework 
and use WKB waves for describing a. Then, in order for the two versions of transition 
amplitudes to coincide, the Wronskians of the WKB waves are univoquely fixed. There¬ 
fore, in the next sections, the gravitational waves shall have unit Wronskians, since we 
must recover the Schrodinger equation in the limit in which a can be described at the 
background field approximation. 


3 The adiabatic treatement in quantum cosmology 

3.1 When matter energy is a conserved quantity 

As a warming up, we first analyse the WDW equation in the simple case when the matter 
energy is a constant of motion (which is a stronger condition than being adiabatically 
conserved). This implies that states with different matter energy are completely decor- 
related. In the next sections we shall drop this assumption and see how the coupling 
to gravity introduces the notion of evolution through interferences among neighbouring 
matter states characterized by adiabatically conserved energy. 

We work in a closed Robertson-Walker universe whose metric is given in eq. (p3|). 
Then, only one constraint survives corresponding to the reparametrisation invariance 
t = t(t'). Classically, this constraint is 

_(~<2 2 _ 2 I A 4 

n = - a — a - + n, a) = 0 (51) 

where a is the scale factor now promoted to a dynamical variable, n a its conjugate 
momentum, Hm is the matter hamiltonian, with 0 and II canonically conjugate matter 
degrees of freedom. A is the cosmological constant. It shall play no other role than to 
participate to the determination of 7i a in terms of a and Hm- 

For the matter energy to be a conserved quantity Hm must have the form 


H m (<4, n, a) = a r h M (4>, II) 


(52) 


(or a sum of such terms). The simple example we consider a very massive field at rest (i.e. 
k — 0). Then H M = h M = (II 2 + M 2 0 2 )/2 where M is the mass of the particles, see ||. 
Here n and 0 are rescalled by a factor a 3 / 2 with respect to the fields defined in section 2.3. 
This rescalling replaces the adiabatically invariant states characterized by instantaneous 
occupation numbers given in eq. (|33|), by truly invariant states characterized by fixed 
occupation number. 

In quantum cosmology, the physical states S(a, 0) are solutions of the WDW con¬ 
straint 


7-fS(a,0)= CA<9 2 — a 2 + Aa 4 + 2GaH M {(j),id<i>, a) S(a, 0) = O 


(53) 


(We neglect any ambiguity due to ordering problems. These would add terms propor¬ 
tional to G 2 /a 2 which could easily be taken into account.) When, Hm is of the form eq. 
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(|52|), Jim is independant of a and commutes with 7i. We can thus work in a basis of 
eigenstates of Iim 


^M\'4^n) E n \'ijjn) 
(t/VnJ'ty'U) = 3m,n 


(54) 


Notice that these eigenstates coincide with the stationary eigenstates of the Schrodinger 
equation id t \i]){t)) = h M \ijj(t)), see section 2.1. Using these states, we can decompose S 
as 

H(a,0) = Y C n^n(a)((f>\i>n} (55) 


Thus, since d a \ip n ) = 0, the different terms in eq. ([55]) are completely uncorrelated. The 
n-th wave function T n (a), entangled to the matter state | ijj n ), obeys 


(~G 2 dl - a 2 + Ao 4 + 2 GaE„) W„(a) = 0 


(56) 


It corresponds to the equation for an oscillator whose time (= a) dependent frequency 
is given by 

a; 2 (a) = (-a 2 + A a 4 + 2 GaE n ) /G (57) 

The solutions of this equation have been discussed in section 2.4. 

To make contact with former works [|3[]|4[] |30| , we recall how a time parameter can be 
extracted from wave packets made out WKB solutions of eq. (|56|). Following^. ||], we 
shall then recall why interactions leading to transitions must be added in order to give 
physical substance to the expressions. 

The WKB solution of eq. (B3) with positive unit wronskian is 


Xn(a) = 


3 —i j a da p n (a ) 

\/ 2 P»(“) 


(58) 


where the momentum p n {a) is the solution of eq. (^) with Hm = E n , hence given by 

p(a, E n ) = p n (a) = —y/—a 2 + A a 4 + 2 GaE n (59) 


These WKB solutions correspond to expanding universes. Their complex conjugate 
describe contracting universes. 

Consider now a superposition, eq. ([55|), made out of WKB solutions with positive 
Wronskian only. Lets further assume that the coefficients C n form a well peaked en- 
veloppe centered around the mean value n such that the phase shifts appearing in the 
sum can be evaluated to first order in n — n only and that the dependance in n — n of 
the prefactors p n (a)” 1 / 2 be neglected, see subsection 3.3 for more details. Then eq. (|55|) 
can be factorized as 


:(a, 0 ) = 


f a dapn(a) 


2 Pfi{o) 


Y, C n e-^n-Efi) d E f a dap(a,E) 


(60) 


16 












The factor of E n — E n 


d E j da p(a, E)\ E = En = tn(a) 


(61) 


is the lapse of proper time evaluated in the classical geometry driven by the mean matter 
energy in virtue of Hamilton-Jacobi equations. 

Thus, well peaked wave packets in quantum cosmology factorize into a gravitational 
piece characterized by the mean matter energy and and a linear superposition of matter 
states multiplied by the conventional phase factors parametrised by the mean time tn{a). 
However, in this model characterized by constants of motion, this quantum superposition 
is physically empty since there is nothing to probe the evolution, i.e. to give physical 
meaning to the phase factors which arise precisely from the comparison of neighbouring 
gravitational waves entangled to the corresponding matter states. Therefore, neither a 
physical interpretation can be given to the coefficients C n nor can a principle be invoqued 
to normalise them. 

On the contrary when the matter hamiltonian no longer commutes with the WDW 
constraint (i.e. which now depends on a) the interactions among matter states gives a 
precise physical meaning to both the phase A E J dadEp(a, E) = A Et and to the coeffi¬ 
cients C n whose normalisation will be univoquely fixed. Indeed, A Et is the phase which 
accumulates between successive transitions and therefore which governs their rate. And 
the coefficients C n obey a Schrodinger equation only if the Wronskians of the correspond¬ 
ing gravitational waves y n (a) are all the same, see the end of section 2.4. 


3.2 Non adiabatic transitions in quantum cosmology 

In this subsection, the matter hamiltonian depends explicitly on a. Examples of such 
Hamiltonians where analysed in section 2.2 and 2.3. They will be further discussed 
in section 4. The coefficients C n then depend on a. To obtain the equation which 
govern their evolution, we must extend the usual adiabatic procedure given in section 
2.1. Indeed, because of the second order character of the WDW equation, a double 
adiabatic development is now required to obtain a linear equation in d a for the C n , First, 
as in section 2.1, we diagonalise instantaneously the matter hamiltonian and secondly, 
we apply an adiabatic decomposition of the type given in eq. ([43]) to the gravitational 
wave functions entangled to the instantaneous matter eigenstates. Moreover, in order 
to determine univoquely the coefficients C n , an additional mathematical condition must 
still be supplied. We proceed formaly by imposing that S(a, (p) satisfies an equation that 
generalyses eq. ([44|). The resulting equation for the C n describes both non-adiabatic 
transitions between matter eigenstates and the backscattering from expanding universes 
into contracting ones. The physical meaning of this equation will be analyzed in the 
next sections. 

Thus, first, we carry out an instantaneous diagonalisation, i.e. at fixed a, of Hm'- 

H M (a)\ip n (a )) = E n {a)\^ n {a)) 

{ip m {a)\ip n (a)) = 8(m, n) (62) 
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We emphasize once more that this diagonalisation does not require the “existence” 
of a Schrodinger equation, although the eigenstates |-0 n (a)) and the eigenvalues E n (a) 
coincide with ones that one would have obtained if one had started, as in section 2.1, 
from the Schrodinger equation whose r.h.s. contains the hamiltonian given by 77 m («(£)). 

Using these instantaneous eigenstates, S(a, 0), solution of the Wheeler-DeWitt equa¬ 
tion (|53|) can be decomposed, as in eq. (|55D but without any lost of generality, as 

(0|S(a)) = E(a,0) = ^^ n (a)(0|'0 n (a)). (63) 

n 

The difference with eq. fl55[) is that the waves <^ n (a) are no longer decorrelated since the 
instantaneous eigenstates |0 n (a)) do now depend on a. In section 3.4, the correlations 
among the ip n (a) will be related to the non-diagonal matrix elements of d a vr^, the kinetical 
term of gravity viewed as an operator acting in the matter Fock states. 

Up to now, no difference exists between the present development and the adiabatic 
treatement of section 2.1. The difference arise from the fact that the WDW constraint 
is second order in d a ■ Therefore, to obtain a first order equation for the coefficients C n , 
we must proceed to a second adiabatic development. To this end, we express each p n (a) 
in terms of the WKB waves Xn(a), eqs. (|58|, |59|) with E n replaced by E n {a ): 

<p n (a) = C n (a ) Xn{a) + V n {a) x* n (a ) (64) 


In the limit where both the adiabatic approximation for the matter states and the WKB 
approximation for gravity are valid, the coefficients C n and T> n are constants. In that 
case, we recover the situation presented in section 3.1 wherein there is no correlations 
among the coefficients C n and T> n . 

To describe these correlations in the general case, we first need to complete the defini¬ 
tion of C n and T> n . We proceed mathematically. Notice however that this mathematical 
choice has deep physical consequences since we shall attribute a physical interpretation 
to the C n . We generalyse eq. (]44|) and require 

(' ^n(a)\id a 15(a)) = p n (a ) [ C n (a ) y n (a) - V n (a) x* n (a)} (65) 

The reason for this choice is that it ensures that the total current carried by S(a, 0) 
contains no terms proportional to d a C n (a ) or d n V n (a). Indeed, one has 


(S(a)|h9 a |S(a)) = f d(j) E*(a, 0)i<9 a S(a, 0) = I C n (a)\- - \T> n (a)\ 2 = C (66) 

** n 


Notice that the absence of relative factors in the above sum directly follows from our 
choice of identical (unit) Wronskians. The same remark also applies to the next two 
equations. We are now in position to determine the a-dependance of C n and V n . Taking 
the derivative of eq. (p3|) and using eq. (|U3|) yields 


daCnXn + 9 a V n X* n ~ C n Xn + E>nX* n ) + (V'n| d a ^ m ) (C m Xm + ^mXm) = 0 (67) 

^Pn m 
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Then taking the derivative of eq. (]65|), inserting it in the WDW equation, and using eq. 
(|67|) to eliminate either d a C n or d a V n yields 


f) C 


Y (dM 

m+n ^y/PnPm 

m 2y/PnPm 

+ daPn e -2i f a Pnd^ 

2 Pn 


( 68 ) 


and the same equation with C n V n , i —i. 

This is our central equation. We emphasize that it is equivalent to the original WDW 
equation (|5lD. Its physical meaning is investigated in the next sections. 


3.3 Time evolution and unitary evolution 

To begin the analysis of eq. (| 68 D , we assume in this section that the coupling between 
expanding and contracting universes may be neglected. In the next section, we shall 
take this unusual quantum coupling into account. 

Thus, we drop the second and third term on the r.h.s. of eq. (|65|). This yields 

d a C n = E e -if a ( Pn -Pm)da ^ (69) 

m/n Z \JPnPm 


In what follows, we shall proceed to the interpretation of this equation in three steps by 
progessively relaxing restrictions imposed on the coefficients C n . 

When comparing eq. ( |69| ) to eq. (0), we see that the main difference is that, in 
eq. (|69|), there are additional dependances in the matter states through the momentum 
of gravity p n {o). Therefore, in a first step, we assume that the C n form a well peaked 
wave packet centered around Cn such that a first order expansion in n — n is a valid 


approximation. Then, using eq. (| 6 lD , we obtain 


d a C n = ( d a^n 

rriy^n 


\<P„) c„ 


(70) 


This equation is identical to eq. More precisely, through t n (a), the proper time 

evaluated in the mean universe, we recover the Schrodinger equation governing non adi¬ 
abatic transitions among instantaneous matter states if one identifies the coefficients C n 
with the probability amplitudes c n to find the matter states at t n (a). We emphasize the 
a posteriori character of this identification. Indeed, it is based on the comparison of eq. 
(0) with the resulting simplified equation governing the C n . We recall that the definition 
of the coefficients C n required three equations, namely eqs. (|63|, [>5[). Had we made 

other choices, we would not have obtained this straitforward identification. Moreover, 
two approximations were also necessary: first that the coupling between expanding and 
contracting universes be neglected and secondly that the C n form a well peaked distri¬ 
bution. For this later condition to be met, we recall that the universe be macroscopic, 
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see H H, The generic character of the recovery of eq. ([7j) from eq. (|68| ) through 
these approximations should be re-emphasized. This whole scenario arises whenever 
one considers how to reach the background field approximation from a more quantized 
framework, see the end of section 2.4 and [11], |I2| . 

As in the cases treated in these references, eq. ( |5P| ) contains higher order terms in 
n — n which define the deviations from the Schrodinger equation when the C n form a well 
defined wave packet. This condition must be met in order to have physically meaningfull 
expansions around mean values. Eq. ([k|) differs from eq. (jTOj) through both the complex 
phase and the prefactors. Both give rise to gravitational corrections since the differences 
arise from the quantum backreaction expressed here by the dependance in n of p n (a). 
To analyse these corrections we define new coefficients C n = C n e~ Pnda which allow to 
rewrite eq. (^9[) as 


^9 a C n ^ [ ( &mnPn T ^ I An) 


Pn+Pr, 
2 V PnPr 


c 




(71) 


where we have defined, as in eq. (||), an effective hamiltonian giving the evolution of C n . 
We now expand the p n (a ) around the mean energy E(a). to obtain 


H-nm ( a ) fimnPn (®) 

+h mn {E n - E)d E Pn + iidai’ml'lpn) 

+8 mn \(E n - Efdlpn + (E n - E m ) 2 

2 °Pn 

+0(AE 3 ) (72) 

The first line determines the overall phase of the C n and encodes no physics. The 
second line contains the linear deviation in energy and gives the Schrodinger equation 
in the mean geometry, eq. The third line contains quadratic terms which furnish 

the gravitational corrections to the energy levels and the instantaneous eigenstates. It 
consists of two terms. The first term arises because the energy of matter modifies the 
propagation of gravity, and hence the phase accumulated by the wavefunction. It can 
be expressed as a shift of the energy levels 

A (K) ~ 2 aS (B " ~ E) (73) 


where in the second equality we have evaluated the correction for a matter dominated 


universe in which p E ~ \JaE/G. (Notice that plays a role similar to the specific 
heat of a large system, see App. B of [||.) The second term of the third line in eq. O, 
is the first gravitational correction to the non adiabatic coupling. By diagonalising H e ff 
one can calculate how it modifies the energy levels. This shift will be much smaller than 
A Grav (E n ) of eq. (|73]) because it carries an additional factor 1/p. 

Because of these corrections, there will be a maximum time t max = 0(E 2 /(E n — E) 2 ) 
after which errors in the phases of the C n (a) are of order 1. Then the evolution of the 
C n (a ) determined by the WDW equation no longer coincides with the one of the c n (t), 


20 









solutions of the Schrodinger equation. The shift in the energy levels will also modify 
the non adiabatic transition amplitudes through the modification of the saddle point 
amplitude. At this point, it is to be emphasized that the calculation of these effects in 
specific cases is best carried out directly through the effective Hamiltonian eq. 0 ) , or 
with eq. (pfj|). Indeed the expansion eq. (| 7^) depends explicitly on the mean energy E 
determining the background, whereas the transition amplitudes are independent of E. 
Therefore quantum cosmology tells us that the (traditional) understanding of physics as 
describing matter processes in a given gravitational background should be considered as a 
convenient approximation with no dynamical justification. This point will be illustrated 
in section 4 by some examples and further clarified in what follows. 

We note indeed that the validity of eq. ([k|) puts no restriction on the C n since 
it requires only the decoupling between expanding and contracting universes. This is 
sufficient to guarantee that J2 n \Cn(a)\ 2 = constant. Therefore one must still interpretate 
C n (a ) as the amplitude of probability to find matter in the n-th state at a, even if they 
do not form a well defined wave packet [| 8 |. In this case of arbitrary distribution of 
C n , it is completely meaningless to insist in using a mean background as a reference 
to parametrize transitions. Thus this raises the interesting question of whether one 
should require at all that the C n form a well peaked wave packet in order to recover the 
notion of a temporal evolution. To this end we note that if one assumes that the matter 
hamiltonian is described by a held theory with a cutoff (at the Planck energy), non 
adiabatic transitions never occurs between states whose energy differs by more than one 
Planck mass. Thus components of 5 with very different matter energies will essentially 
never interact. This offers the possibility that various components of the total waves 
decohere and evolve with their own time. We will not attempt further discussion of this 
issue here. Note however that this problem possesses an analogue in the emergence of 
a temporal description of electronic transitions in atomic collisions in which the atoms 
are described by broad wave packets, see the introduction. 


3.4 Backscattering effects 


We now turn to the coupling between expanding and contracting universes, that is the 
coupling between the coefficients C n and V n in eq. (^). This coupling is inevitable 
because the WDW equation is second order in d a . It also occurs in the simple harmonic 
oscillator problem discussed in section 2.4. However in the present case the situation 
is more complicated because there are backscatterings which keep the matter state un¬ 
changed, i.e. the third term in eq. (^), and backscatterings which change the matter 
state, i.e. the second term in eq. (^). All of these couplings however have the same 
origin: they can all be expressed in terms of matrix elements of 

7 r 2 (a; 4 >, n) = (Aa 4 — a 2 )/ + 2aiAw(a; 0, n) (74) 


viewed as an operator parametrised by a and acting on matter states (here / is the 
identity operator on these states). Indeed, using eq. (H), eq. (|68D can be written as 

1 


ld a C n = 


m^n 


2( Pn Pm) yjPnPi 


-C 
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(75) 


+ J2 i (' l l ; m\daTT 2 a \^n) 

m 


2 ip n + Pm)y/PnPm 


V 


m 


The first term in this equation corresponds to the first term in eq. (^), and the second 
term regroups the two last terms of eq. fl68f) . Therefore, the operator <9 a 7r^ generalyses 
the operator d t H(t) in the Schrodinger case, see eq. (0j), in that it allows this very simple 
rewriting of the WDW equation. 

The coupling between forward and backward propagating universes described by eq. 
(|75|) gives rise to new effects. First it predicts the generation of backward propagating 
waves, i.e. T> n ^ 0, from initialy purely forward waves and the simultaneous increase 
the norm of the “forward” coefficients C n so as to ensure conservation of the Wronskian 
eq. (|66| ) . This is strictly the analogue of the Klein paradox in relativistic quantum field 
theory. Indeed, no positive local)] conserved quantity can be either constructed from 
the solutions of the Wheeler-DeWitt, contrary to the situation in quantum mechanics. 
It is possible that the resolution of this problem requires that one proceeds to third 
quantization wherein S(a) is interpreted as an Heisenberg operator acting on universe 
Fock states. 

Secondly, eq. ([7^) stipulates that there are quantum interferences among matter 
states associated with the expanding and contracting solutions. This means that the 
knowledge of the coefficients V n (a ) is required to determine the evolution of the C n (a) 
at other values of a. Thus the set composed of the C n only form an open system. 

In a future publication we hope to return to these aspects of quantum cosmology. 


4 Applications of the formalism 


In this section we return to the examples of non adiabatic transitions of section 2. We 
shall now treat them in the context of quantum cosmology when gravity is described in 
the adiabatic approximation hence with the amplitudes of matter transitions described 
by eq. . The essential difference is that whereas in section 2 it was merely convenient 
to parametrise the evolution by a, see eqs. ( 1 , 0 ), here it is mandatory to do so. This is 
particularly clear when calculating gravitational corrections to the transition amplitudes, 
since these corrections depend only on the structure of the energy levels E n (a), and 
are independent of the mean energy E and the mean time t defining the reference 
background. 

However before starting we must address the following problem. In the examples 
treated in section 2, the particles have non vanishing momentum which will backreact 
onto gravity and breack the symmetry of minisuperspace. The way around this difficulty 
is to note 10 that if the momentum of the particles is small, then one can treat the 
corresponding metric fluctuations perturbatively. Thus one expands both of the spatial 
metric h tJ and the lapse and shift in a fourier series, and keeps in the Einstein-Hilbert- 


See 117 for recent attemps to build non-local positive definite conserved quantity. 
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matter action only quadratic terms with k ^ 0. The Hamiltonian then becomes 


J d 3 xN^H^ 


N°(k = 0)H o (k = 0) + J2 N^H^-k) 

0 


( 76 ) 


All the new constraints H^(k) which are multiplied by N^^k ^ 0) are linear in the 
gravitational perturbation hjj(k) and simply serve to eliminate the non propagating 
parts of hij. The only non trivial constraint is Ho(k = 0). Thus to this order, one has 
recovered a minisuperspace model, but which contains matter and linearised gravitons 
with non vanishing momenta. Since they are not coupled, one can treat them seperatly. 
We shall for simplicity only consider the matter part of the Hamiltonian. 

We first consider the recovery of the Golden Rule in quantum cosmology. The in¬ 
stantaneous eigenstates | ijj±) and eigenenergies E± where already obtained in section 
2.2. Therefore the effective hamiltonian eq. m for the coefficients C± multiplying \ip±) 
is 


H eff (a ) = 


P+( a ) 

-W- 


i{d a ip_ |'0+)f^= 5 \ / p + (a) i.(d a v^ v .) 

p-(a) J \ -i(d a i/}-\ip+) p-(a ) 

(77) 

where in the second equality we have neglected the gravitational corrections to the non 
adiabatic couplings. Here the momenta of gravity p± characterizing the gravitaional 
waves entangled to the matter states |are 


P± 


^2 Ga(E T + E ± (a)) 

\j 2 Ga(E T + ¥ + m ^ + k / a ± ^( M -m- k/a) 2 /A + g 2 /a^) (78) 


We have chosen to work in a flat matter dominated universe with conserved matter 
energy E r . Since eq. (^) has exaclty the same structure as eq. (f|), we can calculate 
the non adiabatic transition amplitude as in section 2. We recall that the probability to 
make a non adiabatic transition is the probability for the excited atom not to decay. In 
the saddle point method it is given by 


Pr 


no transition 


CL* 

— g—2/m f da [p+(a)-p_ (a)] 


(79) 


where a* is the complex value of a where the momenta of the universe are equal 


p+(o*)-p_(a*)=0 


( 80 ) 


This gravitational resonnance condition leads back to the background held resonnance 
condition E + (a*) — E_(a*) = 0, since p±(a) = p(a,E±), see eq. (^). Thus a* is 
unafected by the dynamical character of gravitys and still given by eq. (0)- 
The integral which yields the transition probability is simply 

f a + 

Ik = Im / da [p + (a) — p-(a)\ (81) 

J ao 
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In order to estimate this integral, we first expand the integrant around the average total 
energy E = E r + {E + + E_)/ 2. The important word in the previous sentence is “total”. 
This is why this expansion is physically meaningfull: the mass-energy of the matter 
content of the universe which is not involved in the process determines the “inertia” of 
gravity. Indeed, we obtain 

4 = Im [ da d E p(a,E)\ E=E [4+(a) - E_(a)\ + 0 
J ao 

= Im£ da ^(a) [E+(a) - E_(a)} + O ) (82) 

Note that this is still not strictly equivalent to the background field approximation since 
the “time” parameter t determined by 


(t) 


d E p(a, E)\ e=e — — 


(83) 


depends both on the initial state of energy E e + M and on the final state of energy 
E t + m + k/a. Thus successive transitions, corresponding to successive k , are described 
by (slightly) different time parameters. A single classical time parameter could be chosen 
to consist of the atom in its excited state of energy M, call it t(a). This is particularly 
appropriate in the present case because we are considering the probability for the atom 


to stay in its excited state. Then eq. fl82[) becomes 


r a * dt 

4 = Im / da—(a) [E + (a) — E_(a)] (1 — 
J ao da 

((E + -E.f\ 


AM 


ao) 


2{E e ~\~ A/) ao 


+ O((a-a 0 ) 2 )) 


+0 


E't 


(84) 


The first term corresponds to eq. O evaluated in the background driven by E e + M. 
The second term is the first gravitational correction to the decay probability. It leads 
to a modification of the decay rate of the order of A M/E E - Note that the correction 
decreases with increasing ao = k/AM. Indeed as the universe expands it is larger and 
larger and the background approximation becomes better and better. 

We now apply the same techniques to the calculation of pair creation amplitudes in 
quantum cosmology. As in section 2.3 we consider only the transition from the vacuum 
state to the state containing one pair of momentum fc, —k. The only relevant quantities 
are the momenta of gravity associated with these two states. They are given by 


Po(ci) = V Aa 4 — a 2 

Pk-k(a) = \j Aa A — a 2 + AGa\Jm 2 + k 2 /a 2 


(85) 


Then, in the saddle point calculation, the amplitude to produce a pair is simply given 
by 


14 


■pair 


2 _ - 2 Im f a da [ Po(a)~Pk,-k(a )] 


( 86 ) 
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The saddle point a* of eq. (|86|) is again determined by the equality of the two mgravi- 
taional momenta: Po(a*) — Pk,-k( a *) — 0. Thus a* is still given by a* = ik/m, which is 
the result obtained in the background field approximation eq. (|38| ). 

To estimate the integral of eq. (^) one can expand the integrant around empty de 
Sitter space, by exploiting the inertial delivered by the cosmological constant. However, 
it is more interesting to note that this expansion around a background is neither phys¬ 
ically intrinsic nor mathematically necessary. Indeed, the integral of eq. ( PB| ) is easily 
evaluated by deforming the integration over a: by first integrating along the real axis 
to a = +oo, then along the circle at infinity from a = +oo to a = +ooe i7r / 2 and fianlly 
along the imaginary axis from a = ioo to a = a* = ik/m. The first piece of the inte¬ 
gral, being real, does not contribute to the pair production probability, the second piece 
yields the background field result, eq. ([36|), and the third piece yields the gravitational 
corrections to the pair production probability. One finds that for large k, corresponding 
to productions occuring far from the classical turning point a = l/\/A, the probability 


is 


\A 


pair | 


~ e 


-+0(Gtra*/fc 3 A 3 / 2 ) 


(87) 


Thus in minisuperspace, the gravitational corrections to pair creation decrease as the 
time (which is related to k by the real part of eq. (^)) of production increases. This 
is natural since at late times the inertia delivered by the cosmological constant is huge, 
and thus, accordingly, the background field approximation better and better, ft should 
however be noted that there probably are local gravitational corrections to pair creation 
which survive even at late times. In particular in de Sitter space these should be related 
to changes of the area of the horizon|33|]. Similar corrections should also obtain for the 


radiative transitions discussed in the previous paragraph. 


5 Conclusion 

In this article we have developped a formalism for analysing the Wheeler-DeWitt equa¬ 
tion based on a double adiabatic developpement. The matter part of the wave function 
was expanded in terms of the adiabatic matter states and the gravitational part in terms 
of the adiabatic (WKB) waves. The equation for the coefficients of this devclloppement 
is the main result of this paper. Because of its generality this formalism is usefull to 
study conceptual problems such as the recovery of a temporal evolution for matter, and 
the identification of the probability amplitude to be in a given matter state. We have 
applied it to two examples: radiative transitions and particle creation, and we have 
compared how these processes are described in the background field approximation and 
in quantum cosmology. 
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